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Abstract This letter reports the results of a search for 
top and bottom squarks from gluino pair production 
in 4.7 fb _1 of pp collisions at y/s = 7 TeV using the 
ATLAS detector at the LHC. The search is performed 
in events with large missing transverse momentum and 
at least three jets identified as originating from a b- 
quark. Exclusion limits are presented for a variety of 
gluino-mediated models with gluino masses up to 1 TeV 
excluded. 

Supersymmetry (SUSY) [THS] provides an extension 
of the Standard Model (SM) which resolves the hi- 
erarchy problem |10H13j by introducing supcrsymmct- 
ric partners of the known bosons and fermions. In the 
framework of the i?-parity conserving minimal super- 
symmetric extension of the SM (MSSM) Q3HH], SUSY 
particles are produced in pairs and the lightest super- 
symmetric particle (LSP) is stable, providing a possible 
candidate for dark matter. In a large variety of models, 
the LSP is the lightest neutralino (xi)- The coloured 
superpartners of quarks and gluons, the squarks (q) 
and gluinos (g), if not too heavy, would be produced in 
strong interaction processes at the Large Hadron Col- 
lider (LHC) and decay via cascades ending with the 
LSP. The undetected LSP results in missing transverse 
momentum - whose magnitude is referred to as E™ lss 
- while the rest of the cascade yields final states with 
multiple jets and possibly leptons. In the MSSM, the 
right-handed and left-handed squarks, q R and q L , can 
mix to form two mass eigenstates q\ and q\ . The mixing 
effect is proportional to the masses of the SM fermion 
partners and can therefore be large for the third gener- 
ation. This may lead to the lightest sbottom (b\) and 
stop (tjj mass eigenstates being much lighter than the 
other squarks. As a consequence, b 1 and t x could be pro- 
duced with relatively large cross-sections at the LHC, 



cither directly in pairs, or through gg production fol- 
lowed by g — > b\b or g — > t x t decays. 

This letter extends the search for gluino-mediated 
bi and t x production at ATLAS reported earlier 119] . 
which used 2.05 fb _1 of data collected in the first half 
of 2011 at a centre-of-mass energy of 7 TeV. The present 
analysis comprises the full 2011 dataset of 4.7 fb _1 and 
adopts an improved selection that requires large E^ lss , 
no electron or muon and at least three jets identified 
as originating from 6-quarks (6-jets) in the final state. 
Results are interpreted in four simplified models where 
sbottoms or stops are the only squarks produced in the 
gluino decays, leading to final states with four &-quarks. 
Searches in similar scenarios have also been reported by 
the CMS Collaboration [20]. 

The ATLAS detector [21] consists of inner track- 
ing devices surrounded by a superconducting solenoid, 
electromagnetic and hadronic calorimeters and a muon 
spectrometer with a toroidal magnetic field. The inner 
detector provides precision tracking of charged particles 
for 1 77 1 < 2.E0. It is immersed in a 2T magnetic field 
from the solenoid and consists of a silicon pixel detec- 
tor, a silicon strip detector and a straw tube tracker 
that also provides transition radiation measurements 
for electron identification. The calorimeter system cov- 
ers the pseudorapidity range |r?| < 4.9. It is composed 
of sampling calorimeters with either liquid argon (LAr) 
or scintillating tiles as the active medium. The muon 

1 ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The x-axis 
points from the IP to the centre of the LHC ring, and the 
y axis points upward. Cylindrical coordinates (r, <f>) are used 
in the transverse plane, <f> being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 6 as rj = — lntan($/2).The distance AR in the 
77 — 4> space is defined as AR = (Arj) 2 + (A0) 2 . 
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spectrometer has separate trigger and high-precision 
tracking chambers which provide muon identification 
and momentum measurement for \t]\ < 2.7. 

Samples of simulated events are used for the descrip- 
tion of the background and to model the SUSY signal. 
The dominant sources of background come from events 
with 6-quarks in the final state. Monte Carlo (MC) sam- 
ples of tt, W/Z and diboson events produced in asso- 
ciation with light- and heavy-flavour jets are generated 
with ALPGEN [22J and the parton distribution function 
(PDF) set CTEQ6L1 [53]. These samples are generated 
with different maximum numbers of additional partons 
at the matrix-element level. Diboson samples are gen- 
erated with up to three additional partons, it+jet and 
Z(—> £ + £~)+]ets (£ — e,fj,,r) samples with up to five ad- 
ditional partons, and W{— > ^)+jets and Z(— > Pz/)+jet 
samples with up to six additional partons. Single top 
quark production is simulated with MCONLO [24J using 
the next-to- leading-order (NLO) PDF set CTEQ6 . 6 [25] . 
The fragmentation and hadronisation for the ALPGEN 
and MCONLO samples are performed with HERWIG [5^1 
[27] . using JIMMY [28] for the underlying event. Sam- 
ples of ti+W, ti+Z and ti+WW events are generated 
with MADGRAPH [25] interfaced to PYTHIA 30 . The sig- 
nal samples are generated using Herwig++ [3T]. The 
MC samples are processed through the ATLAS detector 
simulation [32] based on GEANT4 33 taking into account 
the effect of multiple pp interactions per bunch cross- 
ing. For the comparison with data, all SM background 
cross-sections are normalised to the results of higher- 
order calculations when available, using the same values 
as in Ref. Q]5]. 

Jets are reconstructed from three-dimensional calorime- 
ter energy clusters using the anti- kt jet algorithm |34[ 
135] with a radius parameter of 0.4. The measured jet 
energy is corrected for inhomogeneities and for the non- 
compensating nature of the calorimeter by using pt- 
and 77-dependent correction factors, and additional cor- 
rections for pile- up are applied [33] • Jets are required to 
have pt > 20 GeV, and are reconstructed in the range 
|?7| < 4.9. Events are rejected if they include jets failing 
the quality criteria described in Ref. |36j , or if there is 
any selected jet with \rj\ < 2 for which the scalar sum of 
the transverse momenta of its associated tracks is less 
than 5% of the jet pt- A neural- network-based algo- 
rithm [3TJ is used to identify jets containing a 6-hadron 
decay. This uses as inputs the output weights of differ- 
ent algorithms exploiting the impact parameter of the 
inner detector tracks, the secondary vertex reconstruc- 
tion and the topology of b- and c-hadron decays inside 
the jet. Three operating points are used, corresponding 
to efficiencies of 60%, 70% and 75% for tagging 6-jets in 
a MC sample of tt events. In all cases the tagging rate 



is less than 2% for light-quark and gluon jets, 10% for r 
leptons decaying hadronically and 25% for c-quark jets. 
The 6-jets are identified within the nominal acceptance 
of the inner detector < 2.5) and are required to 
have pt > 30 GeV. To compensate for the differences 
between the 6-tagging efficiency and the mistag rates 
in data and MC simulation, 6-tagging scale factors are 
applied to each jet in the simulations, as described in 
Refs. [37H3H]. 

Electrons are reconstructed from energy clusters in 
the electromagnetic calorimeter matched to a track in 
the inner detector. Electron candidates are required to 
have pt > 20 GeV and |?7| < 2.47 and must satisfy 
the "medium" selection criteria described in Ref. [40] . 
Muons candidates are identified using a match between 
an extrapolated inner detector track and one or more 
track segments in the muon spectrometer, and are re- 
quired to have pt > 10 GeV and | ^7 1 < 2.4. 

To resolve overlaps between reconstructed jets and 
leptons, jets within a distance of AR — 0.2 of an elec- 
tron candidate are rejected. Furthermore, any lepton 
candidate with a distance AR < 0.4 to the closest re- 
maining jet is discarded. Events containing any remain- 
ing electrons and muons are vetoed in the control and 
signal regions defined in Tables 1 and 2. 

The measurement of the missing transverse momen- 
tum two-dimensional vector (and its magnitude £"™ lss ) 
is based on the transverse momenta of all remaining jets 
with I77I < 4.9, all electron and muon candidates and all 
calorimeter clusters not associated to such objects. 

Events are selected using triggers requiring one high 
Pt jet and i?™ lss . Different trigger thresholds were used 
to cope with the increasing luminosity. These triggers 
are fully efficient for this analysis, which requires one jet 
with pt > 130 GeV and E T liss > 160 GeV at the offline 
reconstruction stage. Events must pass basic quality 
criteria to reject detector noise and non-collision back- 
grounds. They are also required to have a reconstructed 
primary vertex associated with five or more tracks with 
Pt > 0.4 GeV; when more than one such vertex is found, 
the vertex with the largest summed p\ of the associ- 
ated tracks is chosen as the hard interaction vertex. 
Events are required to have at least three 6-tagged jets, 
and two jet-multiplicity regions (Nj > 4 and Nj > 6) 
are considered by selecting jets with |?7| < 2.8 and 
Pt > 50 GeV. 

Two variables are calculated from the reconstructed 
objects to further select the events: m s and A</> m j n . 
The effective mass m e ff is defined as the scalar sum of 
the Kj? lss and the pt of all selected jets in a given jet- 
multiplicity region. The A<^ m i n is defined as the min- 
imum azimuthal separation between the selected jets 
and the missing transverse momentum direction. Plac- 
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Common criteria: lepton veto, p^ > 130 GeV, Common criteria: lepton veto, p^ > 130 GeV, 

> 3 fr-jets, E^ iss /m etf > 0.2, A</> min > 0.4 = 2 b-jets, E™ isa /m c{{ > 0.2, A0 min > 0.4, 
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Nj 


Z7>miss 


m etf 


6-tag OP 


Eh 


? iss > 160 GeV, m e ff > 500 GeV 


SR4-L 


> 4j 


>160 GeV 


>500 GeV 


60% 


CR 


Nj 


b-tag OP 


corresponding SR 


SR4-M 


> 4j 


>160 GeV 


>700 GeV 


60% 


CR4-60 


> 4j 


60% 


SR4-L, SR4-M 


SR4-T 


> 4j 


>160 GeV 


>900 GeV 


70% 


CR4-70 


> 43 


70% 


SR4-T 


SR6-L 


> 6j 


>160 GeV 


>700 GeV 


70% 


CR6-70 


> 6j 


70% 


SR6-L 


SR6-T 


> 6j 


>200 GeV 


>900 GeV 


75% 


CR6-75 


> 6j 


75% 


SR6-T 



Table 1 Definition of the five signal regions based on the Table 2 Definition of the four control regions used to esti- 
number of jets (Nj), the E^ 1SB , m e ff requirements and the mate the tt+jets background, 
b-tagging operating point (OP). 



ing the requirements A0 min > 0.4 and E™ lss /m c ff > 
0.2 reduces the amount of multi-jet background, where 
jrjmiss j-ggvj^g f rom mis-reconstructed jets or from neu- 
trinos emitted close to the direction of the jet axis. 

Two sets of signal regions are defined which yield 
good signal sensitivity for the various models and pa- 
rameter values studied here. They are characterised by 
having at least four (SR4) or six (SR6) jet candidates, 
no electron or muon, and are further classified as loose 
(L), medium (M) or tight (T) depending on the E^ iss 
and m c ff thresholds and on the 6-tagging operating point 
The requirements that characterise each signal region 
are summarised in Table [TJ 

The main source of reducible background is the pro- 
duction of ti events in association with additional jets 
followed by the leptonic decay of one W boson, where 
the lepton is not reconstructed or is misidentified as a 
jet (mainly through the hadronic decays of a r lepton). 
This background is estimated by normalising the MC 
event yield in the signal region to the extrapolated event 
yield observed in a tt-dominated control region. Sys- 
tematic uncertainties that are correlated between the 
control and the signal regions largely cancel out in this 
procedure. Additional sources of reducible background 
are single top, tt+W/Z and TV/Z-I- heavy-flavour jets. 
Their contributions are taken from MC simulation and 
account for 10% to 20% of the total background depend- 
ing on the signal region. The irreducible background 
ti+bb is also estimated from MC simulation and ac- 
counts for about 10% of the total background in all sig- 
nal regions. The reducible contribution from multi-jet 
events is estimated with a data-driven method, based 
on a jet response smearing technique [41] . and is found 
to account for less than 5% of the total background in 
all signal regions. 

Four control regions where the tt+jets background 
accounts for more than 80% of the total yield are de- 
fined by applying the same jet requirements and lep- 
ton veto as in the signal regions, but requiring exactly 
two 6-jets instead of three or more. The requirements 
m GfT > 500 GeV and E™ iss > 160 GeV are applied to 
all control regions to make them kinematically similar 



to the signal regions, while reducing the contamination 
from possible SUSY signal events. The definition of the 
control regions is summarised in Table [5J The numbers 
of expected SM events in the four control regions, as 
predicted by the jet smearing technique for multi-jet 
events and by MC simulation for other processes, are 
compared to those observed in data in Table |3J The 
results agree within experimental errors. 



CR 


tt+jets 


others 


SM 


data 


CR4-60 


330 ± 90 


65 ± 25 


395 ± 115 


402 


CR4-70 


490 ± 125 


100 ± 35 


590 ± 160 


515 


CR6-70 


38 ± 11 


7 ± 3 


45 ± 13 


46 


CR6-75 


40 ± 12 


10 ± 4 


50 ± 15 


52 



Table 3 Expected numbers of SM events and observed data 
events in the four ti control regions. The contribution from 
tt+jets events is taken directly from MC simulation. The 
column "others" includes the contributions from single top, 
ti+bb, ti+W I Z and VK/Z+jets processes, also estimated 
from MC simulation, and the contribution from multi-jet 
events which is estimated with the jet smearing technique 
and accounts for less than 5% of the total background. The 
column "SM" shows the total expected background and is the 
sum of the columns "tt+jets" and "others" . The uncertainties 
include all detector-related systematic uncertainties. 

The dominant detector-related systematic effects are 
due to the jet energy scale (JES) and resolution (JER) 
uncertainties, and the uncertainty on the 6-tagging ef- 
ficiency and mistag rates. The JES uncertainty is de- 
rived from a combination of simulations, test beam data 
and in-situ measurements [36], and includes additional 
uncertainties due to the jet flavour and nearby jets. 
Uncertainties on the JER are obtained with an in-situ 
measurement of the jet response asymmetry in di-jct 
events. These uncertainties on jets are propagated to 
the E™ lss measurement, and additional uncertainties 
on E!f lss arising from energy deposits not associated 
with any reconstructed objects are also included. The 
6-tagging uncertainty is evaluated by varying the rj-, pr- 
and flavour-dependent scale factors applied to each jet 
in the simulation within a range that reflects the sys- 
tematic uncertainty on the measured tagging efficiency 
and mistag rates. 
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The systematic uncertainties in the modelling of the 
tf+jets background are assessed as follows: the uncer- 
tainty due to the choice of the MC generator is esti- 
mated by comparing the leading-order ALPGEN gener- 
ator to the MC@NLO generator; the uncertainty due 
to the factorisation and matching scale ambiguities in 
ALPGEN arc estimated by independently varying their 
nominal settings by factors of one half and two; the un- 
certainty due to the finite number of additional partons 
at the matrix-element level is assessed by comparing in- 
clusive ALPGEN samples generated with up to three and 
up to five extra partons. Finally the PDF uncertainties 
are estimated using the MSTW2008NNL0 PDF set. Uncer- 
tainties of 100% are assumed for the multi-jet prediction 
and for the cross-section of it and W/Z production in 
association with a 66 pair. For ti+W / Z production, an 
uncertainty of approximatively 70% has been derived 
from the variations of the factorisation and renormali- 
sation scales and from the PDF uncertainties 42 . 

The tt+jets yield in each signal region is extrapo- 
lated from the measured number of events in the cor- 
responding control region (as per Table [2]) using a fit 
based on the profile likelihood method [43]. Each pair 
of control and signal regions is fitted separately, assum- 
ing no signal events. The free parameter in each fit is 
the tt+jets overall normalisation scale, while the con- 
tributions from subdominant background processes are 
fixed at the expected value. Systematic uncertainties 
are treated as nuisance parameters constrained with a 
Gaussian function and correlations are taken into ac- 
count where appropriate. The results of the fits and 
the numbers of observed events for each signal region 
are summarised in Table |U The fitted values of the nor- 
malisation factors for tt+jets are compatible with one 
and the main impact of the data-driven estimate is a 
reduction in the uncertainty by approximately a fac- 
tor of two. Figure [T] shows the measured m c g distribu- 
tions and the MC predictions for the SM backgrounds 
in each signal region. Also shown are the prediction of 
two benchmark signal models described below. 

The reliability of the MC extrapolation of the tt 
background to larger 6-jet multiplicities has been checked 
in validation regions defined with kinematic cuts similar 
to those used in the control and signal regions, except 
that exactly one isolated electron or muon is required. 
The transverse mass of the lepton and the E™ ss is re- 
quired to be less than 100 GeV in all validation regions 
to minimise the possible contamination from stop pro- 
duction. The extrapolated event yield in the validation 
regions with at least three 6-jets from the validation re- 
gions with exactly two 6-jets is found to be consistent 
with the number of observed events for all 6-tagging 
operating points. 
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tt+jets 
(MC) 


others 


SM 


data 


SR4-L 


33.3 ± 7.9 
(32.6 ± 15.4) 


11.1 ± 4.9 


44.4 ± 10.0 


45 


SR4-M 


16.4 ± 4.1 


6.6 ± 2.9 


23.0 ± 5.4 


14 


(16.1 ± 8.4) 








SR4-T 


9.6 ± 2.1 
(11.4 ± 5.4) 


3.7 ± 1.6 


13.3 ± 2.6 


10 


SR6-L 


10.3 ± 3.3 
(10.0 ± 6.2) 


2.4 ± 1.4 


12.7 ± 3.6 


12 


SR6-T 


8.3 ± 2.4 
(7.9 ± 5.3) 


1.6 ± 1.1 


9.9 ± 2.6 
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Table 4 Comparison between the results of the fits and the 
numbers of observed events in the five signal regions. The 
tt+jets event yield predicted by the MC simulation is quoted 
in parentheses. The column "others" includes the contribu- 
tions from single top, ti+bb, ti+W /Z, VK/Z+jets and multi- 
jet processes. Multi-jet events contribute less than 5% of the 
total background. The column "SM" shows the total expected 
background and is the sum of the columns "tt+jets" and "oth- 
ers" . The uncertainties include the statistical and systematic 
uncertainties. 



The background predictions have been further vali- 
dated using a data-driven method that simultaneously 
estimates all SM background contributions with at least 
one misidcntified 6-jet. This method consists of predict- 
ing the number of jets originating from 6-quarks in each 
event by solving a system of equations based on the 
number of 6-tagged and non 6-tagged jets in the event, 
along with the 6-tagging efficiency and mistag rates. 
Consistent background predictions with respect to the 
fit results have been found in all signal regions. 

Limits for non-SM signal at 95% confidence level 
(CL) are derived by testing the signal plus background 
hypothesis in each signal region with the CL S prescrip- 
tion [43]. These limits are obtained with fits similar to 
those used to estimate the background in each signal re- 
gion, except that the number of observed events in the 
signal region is added as an input to the fit and a sec- 
ond free parameter for the non-SM signal strength, con- 
strained to be non-negative, is adjusted in the likelihood 
maximisation. This additional free parameter ensures a 
proper treatment of the expected signal contamination 
in the control regions when the results are interpreted 
in the framework of specific SUSY scenarios. Model- 
independent upper limits at 95% CL on the number of 
signal events and on the visible cross-section (defined as 
the cross-section times kinematic acceptance times ex- 
perimental efficiency) for non-SM contributions derived 
for each signal region are given in Table [5] 

These data have been used to derive limits in the 
parameter space of the following SUSY models: 

Gluino-sbottom model: MSSM scenarios where 
the b 1 is the lightest squark, all other squarks are heav- 
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m e „ [GeV] m e „ [GeV] 

Fig. 1 Distribution of m e ff for SR4-L and SR4-T (top) and SR6-L and SR6-T (bottom). The hatched band shows the 
systematic uncertainty on the MC prediction, which includes both experimental uncertainties (among which JES and b- 
tagging uncertainties are dominant) and theoretical uncertainties on the background normalisation and shape. The label 
"others" includes the contributions from single top, ti+bb, ti+W /Z, W / Z+jets and multi-jet processes. The lower plot in each 
figure shows the ratio of the observed distribution to that expected for the SM background. Two signal points (with small and 
large mass splitting between the gluino and the LSP) for the Gbb and Gtt models described in the text are overlaid. 





Obs (exp) 95% CL upper limit 


SR 


N events 


crvis(fb) 


SR4-L 


23.8 (23.4) 


5.1 (5.0) 


SR4-M 


8.6 (12.8) 


1.8 (2.7) 


SR4-T 


7.1 (9.2) 


1.5 (2.0) 


SR6-L 


9.6 (10.1) 


2.0 (2.1) 


SR6-T 


7.1 (8.3) 


1.5 (1.8) 



Table 5 Observed (expected) 95% CL upper limits on the 
non-SM contributions to all signal regions. Limits are given 
on numbers of events and in terms of visible cross-sections 
defined by cross-section times kinematic acceptance times 
experimental efficiency. Systematic uncertainties on the SM 
background estimation are included in the limits. 



branching ratio for g — > b x b decays is 100%. Sbottoms 
are produced via gg or by b 1 b 1 direct pair production 
and are assumed to decay exclusively via b 1 — » b\i, 
where m~ is set to 60 GeV. Exclusion limits are pre- 
sented in the (m~,m- ) plane. 

Gbb model: Simplified scenarios, where b 1 is the 
lightest squark but < m? . Pair production of gluinos 
is the only process taken into account since the masses 
of all other sparticles apart from the x? are se t above 



the TeV scale. A three-body decay via an off-shell sbot- 
tom is assumed for the gluino, yielding a 100% BR for 
the decay g —> bbXv The sbottom mass has no impact 
on the kinematics of the decay and the exclusion limits 
are presented in the (mg,m^o) plane. 

Gluino-stop model: MSSM scenarios where the 
ti is the lightest squark, all other squarks are heavier 
than the gluino, and m 5 > m- + m t , so the branching 

ratio for g — > t^t decays is 100%. Stops are produced 
via gg and t 1 t 1 and are assumed to decay exclusively 
via ^ — > bx% ■ The neutralino mass is set to 60 GeV, 
the chargino mass to 120 GeV and the latter is assumed 
to decay through a virtual W boson. Exclusion limits 
are presented in the (jn~ g , m~ ) plane. 

Gtt model: Simplified scenarios, where t-y is the 
lightest squark but m= < m~ . Pair production of gluinos 
is the only process taken into account since the mass 
of all other sparticles apart from the \i are above the 
TeV scale. A three-body decay via off-shell stop is as- 
sumed for the gluino, yielding a 100% BR for the decay 
g — s- ttxi- The stop mass has no impact on the kinemat- 
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Fig. 2 Exclusion limits in the (m-,m- ) plane for the gluino-sbottom model (top left), in the (m-,m- ) plane for the 
gluino-stop model (top right) and in the (m,g,m^a) plane for the Gbb (bottom left) and Gtt (bottom right) models. The 
dashed black and solid bold red lines show the 95% CL expected and observed limits respectively, including all uncertainties 
except the theoretical signal cross-section uncertainty. The shaded (yellow) bands around the expected limits show the impact 
of the experimental uncertainties while the dotted red lines show the impact on the observed limit of the variation of the 
nominal signal cross-section by lcr theoretical uncertainty. Also shown for reference are the previous CDF | 44II45 |. DO | 46 | and 
ATLAS |T91H51|4T1HH| analyses. 



ics of the decay and the exclusion limits are presented 
in the (m^,m-o) plane. 

The SR4 regions are mostly sensitive to the SUSY 
models where sbottom production dominates, whilst 
the SR6 regions are used to set exclusion limits in mod- 
els characterised by on-shell or off-shell stop produc- 
tion, where top-enriched final states are expected. The 
signal region with the best expected sensitivity at each 
point in the parameter space is used to derive the limits 
at 95% CL. Signal cross-sections are calculated to next- 
to-leading order in the strong coupling constant, includ- 
ing the resummation of soft gluon emission at next-to- 
leading-logarithmic accuracy (NLO+NLL) The 
nominal cross-section and the uncertainty ^Theory are 
taken from an envelope of cross-section predictions us- 
ing different PDF sets and factorisation and renormal- 
isation scales, as described in Rcf. [54 . All detector- 
related systematic uncertainties are treated as fully cor- 
related between signal and backgrounds. In the Gbb 
scenario, the impact of initial-state radiation (ISR) is 
expected to be large in the region with low nig — m^o 



due to the small signal acceptance. Therefore, an uncer- 
tainty on the modelling of ISR is assessed by comparing 
the signal acceptance obtained with the Herwig++ sam- 
ples to the one obtained with dedicated MADGRAPH sam- 
ples generated with additional jets. This uncertainty 
varies from 4% to 35% as a function of m~ g — m^o and 
is included in the ±lcxhoory band. 

The expected and observed 95% CL exclusion lim- 
its in the four models considered above are shown in 
Figure O In the gluino-sbottom model, gluino masses 
below 1000 GeV are excluded for sbottom masses up to 
about 870 GeV using the most conservative — lcxheory 
hypothesis. This extends by approximatively 100 GeV 
the limits derived in the same scenario by the previ- 
ous ATLAS analysis performed with 2 fb -1 [TH] and is 
complementary to the ATLAS search for direct sbot- 
tom pair production, also carried out with 2 fb -1 [47] . 
The exclusion is less stringent in the region with low 
nig — m~ h , where softer jets are expected. Because of 
the kinematic cuts applied, the limits depend on the 
neutralino mass assumption for low mass splitting be- 
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tween the sbottom and the neutralino as shown for the 
Gbb model where gluino masses below 1020 GeV are 
excluded for neutralino masses up to about 400 GeV, 
improving the previous ATLAS limits [19] by approx- 
imatively 100 GeV. In the gluino-stop model, gluino 
masses below 820 GeV are excluded for stop masses 
up to 640 GeV, extending the previous ATLAS lim- 
its [THH1] by approximatively 150 GeV. In the Gtt 
model, gluino masses below 940 GeV are excluded for 
m-r.o < 50 GeV while neutralino masses below 320 GeV 
are excluded for m g = 800 GeV. This search extends 
the exclusion limits on the gluino mass from the AT- 
LAS multi-jet analysis carried out with the same data 
set [48] and from the CMS same-sign dilepton analysis 
performed with 5 fb _1 [3D] by approximatively 60 GeV 
and 130 GeV, respectively, for neutralino masses below 
100 GeV. In the region with low m g — m^o, the limits 
obtained with the CMS analysis are most stringent due 
to the softer kinematic cuts. 

In summary, this letter presents results from a search 
for top and bottom squarks in the decay of gluino pairs 
produced in pp collisions at yfs = 7 TeV, based on 
4.7 ftr 1 of ATLAS data. The events are selected with 
large E™ lss , four or six jets and at least three jets orig- 
inating from 6-quarks in the final state. The results are 
in agreement with the SM background prediction and 
translate into 95% CL upper limits on excluded masses 
for a variety of SUSY benchmark scenarios. Gluino masses 
up to 1 TeV are excluded, depending on the model, 
which significantly extends the previous results. 
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